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Abstract 

A simple pseudo-Hamiltonian formulation is proposed for the linear inhomogeneous sys- 
tems of ODEs. In contrast to the usual Hamiltonian mechanics, our approach is based on 
the use of a non-stationary Poisson bracket, i.e. the corresponding Poisson tensor is allowed 
to explicitly depend on time. Starting from this pseudo-Hamiltonian formulation we develop 
a consistent deformation quantization procedure involving a non-stationary star-product *t 
and an "extended" operator of time derivative D t = dt + • • • , differentiating the *f-product. 
As in the usual case, the *f-algebra of physical observables is shown to admit an essentially 
unique (time dependent) trace functional Tr^. Using these ingredients we construct a com- 
plete and fully consistent quantum-mechanical description for any linear dynamical system 
with or without dissipation. The general quantization method is exemplified by the models 
of damped oscillator and radiating point charge. 

1 Introduction 

The problem of quantum-mechanical description of dissipative systems remains for decades a 
recurrently discussed physical topic with a number of important applications (see PQ-[H| and 
references therein). It has also some theoretical importance as a touchstone for testing various 
quantization methods. Our interest to the problem is inspired by recent developments in de- 
formation quantization [IB]- In this paper we consider the problem of deformation quantization 
for dissipative systems, whose classical dynamics is specified by linear inhomogeneous systems of 
ODEs. 

In spite of a large number of papers devoted to the quantum-mechanical treatment of various 
dissipative systems there is no commonly accepted definition of the dissipation phenomenon itself. 
It seems that the most characteristic property of the dissipation, shared by all the systems known 
as dissipative, is the presence of attractors jTHj. Recall that an attractor A C M is a compact, 
measure-zero subset in the phase space of the system M, possessing the property of being limiting 
set for any trajectory passing through a sufficiently small neighborhood of A. In a simple situation 
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the attractor is just a stable fixed point (sink) or a closed curve (limit cycle); in the general case, 
however, much more complicated dissipative structures may occur (e.g. strange attractors). 

From the Liouville theorem about the conservation of phase-space volume (201 ^ follows im- 
mediately that no smooth compatible symplectic structure can exist in the vicinity of an attractor: 
As the measure of attractor is zero, the Liouville volume form, given by the Pfaffian of the sym- 
plectic 2-form, becomes necessarily infinite at the points of the attractor. The attractors are thus 
obstructions for constructing the conventional Hamiltonian description of the dissipative systems 
and this is also the reason why the terms "dissipative" and "non-Hamiltonian" are often used as 
synonyms. 

In principle, the requirement of smoothness of the symplectic structure is not so crucial for 
deformation quantization as only the Poisson bracket and a trace density are actually used. The 
problem, however, is that even for a linear dissipative system with one-point attractor (e.g. focus) 
the corresponding Poisson bracket appears to be highly nonlinear, making impossible any practical 
calculations [21] (see also [2Zj for similar discussion of Kepler's problem). 

To get round the "no-go theorem" above and obtain a practicable quantization scheme we 
allow the symplectic structure to depend on time explicitly. The idea is as follows: Since the 
phase-space trajectories reach the attractor A only asymptotically (as t — > oo) one can try to 
construct a time dependent symplectic structure Q(t) which would be a smooth function of time 
and such that lim^oo Q(t) = oo. We show that such a non-stationary symplectic structure 
does exist for any linear dynamical system and can be explicitly constructed by the fundamental 
matrix of the linear system. Moreover, for one-dimensional dynamical systems, described by a 
linear second-order equation, one can always find a time-independent Q, whereas for the general 
multi-dimensional system such choice is impossible. In this respect the one-dimensional dissipative 
systems, like the damped linear oscillator, are not indicative examples. 

The unavoidable time dependence of multi-dimensional symplectic structure has further con- 
sequences for nearly Hamiltonian systems, i.e. the systems whose equations of motion have the 
form of a small perturbation over the Hamiltonian ones. No matter how small the value of per- 
turbation is, the perturbed equations may no longer be Hamiltonian, that indicates presence of 
genuine dissipation. It is still possible, however, to construct the first-order action functional if 
one admits a non-stationary symplectic structure; in so doing, the explicit time dependence of Q 
can violate any phase-space polarization [22] (e.g., this may admit no separation of the variables 
onto position coordinates and conjugate momenta) making impossible the passage to the second- 
order Lagrangian formalism. The absence of the second-order Lagrangian description for some 
classical systems with two or more degrees of freedom was recognized long ago [22], [21] , |25| . 
For a nearly Hamiltonian system it is reasonable to require the Lagrangian, if any, to pass into 
the non-perturbed one as the parameter of dissipation (e.g. friction constant) vanishes. The last 
requirement imposes very strong restriction on physically reasonable Lagrangians. For example, 
in [21] we have described all the second-order Lagrangians for the reduced Lorentz-Dirac equation 
in the case of homogeneous magnetic field, but none of them reproduces the free particle motion 
as the charge of the particle tends to zero. The last fact shows clearly the advantage of the 
first-order formalism over the second-order one. For the general discussion on peculiarities of the 
inverse problem of variational calculus in the first-order formalism we refer the reader to [2E] • 

The absence of a phase-space polarization favors also the use of deformation quantization 
over the canonical quantization procedure. Being the function of time, but not of the phase- 
space coordinates, the corresponding Poisson bracket can easily be quantized by the usual Weyl- 
Moyal formula giving rise to a non- stationary star-product * t . It is the point where our approach 
deviates from the conventional scheme of deformation quantization [21], [21] • Notice that in the 
non-stationary case the usual time derivative does not differentiate the * t -product of quantum 
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observables; instead one can define an extended time derivative D t = dt + • ■ ■ that would be 
compatible with the * r product in the sense of the Leibnitz rule. As in the usual case, the * t - 
algebra of quantum observables is shown to admit an essentially unique (non-stationary) trace 
functional Tr t . Using these ingredients, we define the quantum Liouville equation governing the 
evolution of quantum-mechanical states, as well as the rule for computing the expectation values 
of physical observables. As a result we get a complete quantum-mechanical description for any 
linear dynamical system. Moreover, the physical content of the theory is shown to be independent 
of any ambiguities concerning the choice of the non-stationary symplectic structure Q(t) as the 
quantum Liouville equation coincides precisely with the classical one (see Proposition 1). Notice 
that a lot of objections have been published for many years against possibility to consistently 
quantize such systems in the Hamiltonian framework, even in principle, see e.g. |14j-|17j. 

By way of illustration we consider the quantization problem for the damped linear oscillator 
and for the radiating point charge moving in a homogeneous magnetic field. In both cases the 
quantum dynamics seems very reasonable. In particular, the time evolution of mean energy, 
defined in terms of the corresponding unperturbed system, is shown to coincide with evolution of 
the classical energy. 



2 Pseudo-Hamiltonian formulation of linear dynamical sys- 
tems 

We start with an inhomogeneous linear system of ODEs 

x' = A)(t)x j + J\t) (1) 

defined on a linear phase space with coordinates x l . Hereafter the overdot stands for the derivative 
with respect to time t. When A and J are independent of time, one speaks of an autonomous 
system of ODEs. In this paper we are interested in linear dynamical systems which are a (small) 
perturbation of a Hamiltonian one. Although the perturbed system may no longer be a Hamil- 
tonian one 1 , the number of degrees of freedom remains the same. In particular, we will always 
assume the phase space of the system (JIJ) to be even dimensional, i.e. % = 1, ...,2n. 

Our first observation is that any such system can be derived from the variation principle for a 
quadratic action functional if the explicit time dependence is admitted in the integrand. Consider 
the following ansatz: 

S[x] = ^J dtix^ijit)^ - x'Bijtyxf - 2C i {t)x i ) , (2) 

where 

Qij = —Vtji , = Bji , det(fiij) 7^ . (3) 

Structurally, the functional S is similar to the first-order action associated with the Hamiltonian 

H=^x i B ij (t)x i -C i (t)a?, (4) 

but unlike the usual Hamiltonian formalism we allow the symplectic form Q to depend on time. 
This assumption appears to be crucial for the description of multi-dimensional dissipative systems 
as will be seen below. 

this case, the perturbation can't be induced by a perturbation of the Hamiltonian. 
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Varying this action functional, we come to the following equations 
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.^ = x = fr 1 ( b - -ti ] x + n _1 c . (5) 

5x l V 2 / 

In order for these equations to be equivalent to the original ones ((TJ we must set 

A = IT 1 (b - ~f2j , J = fl -1 C, (6) 

or, what is the same, 

ifi=S- OA, C = OJ. (7) 
Decomposing the first matrix equation onto symmetric and anti-symmetric parts, we finally get 

Q = -(nA + A t Q), B=^fiA-A*fiJ, C = OJ, (8) 

A* being the transpose of the matrix A. Only the first relation is nontrivial (it is a linear ODE 
on fi), while the other two are just definitions of the matrices B and C . 

Recall that the square matrix T(t) is called the fundamental solution to the system ((TJ if 

t = AT, T(0) = 1. (9) 

The columns of this matrix constitute the basis in the linear space of solutions to Eqs.Q. Given 
the matrix T, the general solution to the first equation (jSj) can be written as 

n = A t tt A, (10) 

where A = T" 1 , and Qq = — Qq is a constant non-degenerate matrix. The matrix Qq encodes all 
the ambiguity in the definition of quadratic action functional (J2J for the given system of ODEs 

©• 

In the autonomous case, the system ((TJ is explicitly integrable in elementary functions and 
hence the action (j2J can be written in a closed form. By way of illustration let us consider two 
physical examples: the damped linear oscillator, and the non-relativistic motion of a point charge 
with a due regard for the radiation back reaction. 

Example 1. The equation describing the linear oscillator with friction reads 

x + 2ax + u 2 x = 0. (11) 

Here u is the frequency and a > is the coefficient of friction. Introducing the auxiliary variable 

x + ax 
P ~ y/l - a 2 /Z? ' 

one can replace (|llj) with the following pair of first-order equations: 



x = pyl — a 2 juj 2 — ax , p = —u 2 x\Jl — a 2 juS 1 — ap . (12) 



2 Here we use the matrix notation. 
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According to (J2J) and (jHJ), the action functional to this system is given by 



S[x,p] = cjdt(px- 1 - ^/T^J^(p 2 + uA 2 ) + apx'j e 2 " 



(13) 



c being an arbitrary constant. In the regime of aperiodic damping (a > u) the action functional 
becomes complex. 

In this simple case it is also possible to construct a first-order action functional involving the 
canonical symplectic structure. For example, varying the action 

S[x,p] = j dt (px - l -{e~ at v 2 + w 2 e°V)J , (14) 

one gets the Hamiltonian equations 

x = e~ at p , p = -u 2 e at x , (15) 

which are obviously equivalent to Eq.flTTJ). Notice that, contrary to the system (|T2j) . the stationary 
point x = p = of (fTB]) is not an attractor. 

It is significant that both ()13|) and (fT4*j) come to the standard Hamiltonian action for the 
harmonic oscillator when a — > 0. Since the action functional carries all the information about 
classical and quantum dynamics, the last fact makes possible a consistent interpretation of friction 
as a small perturbation over the given Hamiltonian system, rather than something leading to a 
completely different physical system. 

Example 2. The effective dynamics of a non-relativistic charged particle is governed by the 
Lorentz equation [29j 

e 2e 2 

mx = eE + -[x,H] + — x . (16) 

C 6C 6 

Here x(t) 6 R 3 is a trajectory of the particle, E and H are 3-vectors of electric and magnetic 
fields, and the constants c and e denote the light velocity and electric charge of the particle. 
As is seen, the Lorentz equation involves third time derivative of the trajectory. It is the term 
which describes the back reaction of the radiation emitted by the accelerating charge. Since the 
order of the equation is greater than two, it cannot be assigned with a straightforward mechanical 
interpretation: in the realm of Newtonian mechanics, a trajectory of a scalar particle is uniquely 
specified by initial position and velocity. Besides, together with physically meaningful solutions 
Eq.()16|) allows a set of nonphysical ones [2H|. It turns out that both mentioned problems can be 
resolved by means of the reduction of order procedure (see e.g. [21], [2H], EH])- Namely, Eq. ()16l) is 
replaced by a second-order equation x = g(x, x, e) such that all the solutions to the latter would 
solve the former. The last requirement leads to a partial differential equation on the function 
g(x, v, e) having a unique solution with g(x, v, 0) = 0. 

Consider for example a non-relativistic particle moving in a homogeneous magnetic field, i.e. 
E = 0, H = (0, 0, H) and H = const. The reduced Lorentz equation has the form [21] 



x = Ax — By , y = Bx + Ay , z = 



where x = (x, y, z) and 



6 - V6V3 + V9 + 64e 6 F 2 2 4 rr2 „ eH^> 



A= " v " v " - V ' ^~e 4 H\ B — 7 v ^_«eg. (17) 

3 V 3 + V9 + 64e 6 # 2 



ie 2 
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Here we have set m = c = 1. Since the evolution along z represents the free motion and decouples 
from the dynamics in the xy-pl&ne, we restrict our consideration to the first two equations. 

Setting formally A = 0, we arrive at the usual Lorentz equations describing the motion of a 
charge in response to the "effective" magnetic field B = (0,0, B/e). In this case, the trajectories 
are concentric circles. For A ^ the particle spirals at the origin of iq/-plane. So, it is natural to 
regard A as the coefficient of friction. 

In order to construct an action functional for the second-order Eqs. (|17|) , we replace them with 
an equivalent system of first-order ones. Let us choose the auxiliary variables as 

B B 
p = x + —y, q = y-— x . (18) 



Then 



B B 

2 V > y = q+ r 



B B 2 7 B \ B B 2 ( B 

P = ~^Q- + Alp- —y 1 , q = — p- — y + A I q+ —x 



(19) 



Applying the general formulas (J2J) and (JBJ), we arrive at the following expression for the first-order 
action functional: 

S[x, y, q,p] = 4^2 + ^2) / dte' At [2a(t)(px - xp + qy - yq) (20) 
+2b(t)(qx — xq + yp — py) + 2c(t)(pq — qp) + 2d(t)(xy — yx) 



+e{t){p 2 + q 2 ) + f{t){x 2 + y 2 ) + g(t)(px + qy)+]{t){qx - py)} , 

where 

a(t) = A 2 cos(Bt) + ^B 2 (e- At + e At ) , b(t) = A 2 sm(Bt) - ABe At + AB cos(Bt) , 

c(t) = e~ At B + 2A sm(Bt) , d(t) = ^B 3 ( e - At - e At ) - ^B 2 A sm(Bt) + A 2 B(cos(Bt) - e At ) , 

e(t) = e~ At B 2 + A 2 cos(Bt) + A sin(Bt)B , 
f(t) = ^B (B 3 e- At + BA 2 cos(Bt) - A sm(Bt)[B 2 + 2A 2 }) , 

g{t) = -A cos(Bt)B 2 - A 3 cos(Bt) , j(t) = -A 3 sin(Bt) + e~ At B 3 + A 2 B cos(Bt) . 

In the limit of zero friction A — > 0, this action tends to the usual action for the charged particle 
in a homogeneous magnetic field B. 

Notice that the symplectic structure entering the Lagrangian (J2*U|) does not possess x?/-polarization 
(i.e. the 2-form Q does not vanish upon restriction on the 2-plane x = const, y = const) when 
A 7^ 0. This makes impossible the algebraic elimination of p and q from the action pljl and thus 
obtaining a second-order action in terms of x and y. The latter fact agrees well with the general 
statement of Ref.j2I| about nonexistence of a second-order action functional for Eq. ()17j) . which 
would pass to the standard action functional for a free particle when e — > 0. 

It can also be shown that unlike the previous example, the system (|17|) does not admit a 
first-order action involving a stationary symplectic structure and having the standard free limit. 
The unavoidable time dependence of the symplectic structure may be thus viewed as a specific 
feature of mult i- dimensional dissipative systems. 
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Using the Hamiltonian (j3J) and the symplectic 2-form Q it is possible to rewrite Eqs.Q in a 
pseudo-Hamiltonian form. To this end we introduce the following nonstationary Poisson brackets: 



{ F, G}( = n«( t )g-g-, n = n-, (21) 

where F and G are functions of the phase-space coordinates x l and the time t. Clearly, these 
brackets satisfy all the properties of the Poisson brackets: bi-linearity, skew-symmetry, the Leibnitz 
and Jacobi identities. In accordance with (jSJ the evolution of an arbitrary physical observable 
F(t, x) is described by the following equation: 

dF 

— = D t F + {F,H} t , (22) 



where 



H = -x'-ByZ* + C,(t)x l (23) 



and 



is the "extended" partial derivative in t. When f2 is constant, the last term in (jzljl vanishes and 
we arrive at the conventional Hamiltonian equations w.r.t. the canonical Poisson brackets and the 
Hamiltonian ()23j) . The main property of the extended time derivative D t is that it differentiates 
the non- stationary Poisson brackets (|2ip. i.e. 

D t {F, G} t = {D t F, G} t + {F, D t G} t , (25) 

for any F(t, x) and G(t, x). Using this property one can deduce the Poisson theorem: 

dh dG d , „ „, 

(the Poisson bracket of two conserved quantities is conserved). 

In the full analogy with the conventional Hamiltonian mechanics, a state of the system is 
described by a classical distribution function p c i(x,t) subject to the normalization condition 

d/j,p c i(x,t) = 1 , (27) 



dfi = v det Q d 2n x being the Liouville measure associated to the non-stationary symplectic form 
Q. The time dependence of p c i(x,t) is determined by the modified Liouville equation (cf. Eq. 



D tPd = {H, Pd } t . (28) 

Note that the time evolution preserves the normalization condition (|27jl . Indeed, using the obvious 
identity 

r dp{F,G} t = 0, (29) 



where one of the functions F(x) and G(x) vanishes on infinity, one can find 

j t J dpp d (x, t) = - J dp Pcl {x, t) = 0, A = VdetQ. (30) 
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The pure states of the classical system correspond to the ^-distributions 

Pd (x,t) = A- 1 (t)5 2n (x-x (t)), (31) 

supported on the integral trajectories Xo(t) of the system (JTJ). 

As the final remark, let us note that the above pseudo-Hamiltonian formalism is applicable 
for arbitrary (not necessary quadratic) Hamiltonians and can easily be derived/justified in the 
formalism of constrained dynamics (31] applied to the first-order action (J2J). 

3 Deformation quantization of pseudo-Hamiltonian sys- 
tems 

Example 2 of Sec.l shows that the most characteristic feature of multi-dimensional dissipative 
systems is the lack of a phase-space polarization compatible with dynamics: Though at each time 
moment to one can split the phase-space variables x % = {q a ,Pb) on "coordinates" and "momenta" 
satisfying the canonical Poisson bracket relations 

{g a , q b U = , {p«, Pb }t = o , {<f , PbU = S a b , (32) 

these relations may not hold true at the next time moment. This is due to the explicit (and in 
many interesting cases unavoidable) time dependence of the symplectic form Q entering the first- 
order action (|2*|). The absence of a natural polarization favors the use of deformation quantization 
over the canonical quantization procedure. 

In the approach of deformation quantization the classical observables (i.e. functions in the 
phase-space variables) are identified with symbols of operators |2Zj, |2Hj; in so doing, the pointwise 
multiplication of functions is replaced by an associative noncommutative star-product. Using the 
non-stationary Poisson bracket (|2*T|). we define the time dependent * t -product by the usual Weyl- 
Moyal formula 

(F * t G)(x, t) = exp ^IP»"( t )_^_^ F (t, x)G(t, y)\ x=y = F-G + ^{F, G} t + 0(h 2 ) , (33) 

It is clear that the time dependence of II does not affect the associativity of the Weyl-Moyal 
star-product, so we have 

(F * t G) * t H = F * t (G * t H) , VF,G,H e C°°(M 2n ) . (34) 

In order to define the notion of a quantum state we endow the * f -algebra with the following 
trace functional: 

Tr t (F) = — J d^F{x) . (35) 

The basic property of the trace (specifying it up to multiplication on an arbitrary function of t) 
is vanishing on * r commutators, i.e. 

Tr t ([F,G\ t ) = 0, (36) 

where [F,G] t = F * t G — G * t F and at least one of the functions F and G vanishes on the 
infinity together with all its derivatives. Actually, the identity (J36)) follows from the stronger one: 
Tr t (F * t G) = Tr t (F ■ G). 
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A pure state of a quantum mechanical system is described by a Wigner function p(t, x) subject 
to the following conditions: 

p* t p = p, Tr t (p)=l. (37) 

The quantum counterpart of the classical Liouville equation, governing the evalution of a quantum 
state, reads 

ihD t p + [p, H] t = 0, (38) 
where the extended time derivative D t , defined by Eq. (JMJ), can also be written as 



dF 1 • 

D t F = — - ^ij(x l * t [x\ F\ t +[x>,F ] t * t x l ) . (39) 

Note that Eq. (J38|) makes sense for an arbitrary (not necessary quadratic) Hamiltonians. For a 
constant Poisson bracket, (|38|) reproduces the usual von Neumann's equation for the symbol of 
statistical operator p. A simple direct calculation shows that the operator D t differentiates the 
* r product: 

D t {F * t G) = (D t F) * t G + F* t (D t G) . (40) 

As a consequence, every solution to Eq. (|58j). satisfying the idempotency condition (}3*T|) at some 
initial time moment, will satisfy this condition in all subsequent time moments. The specific choice 
of the integration measure in the definition of trace functional (|35|) provides the conservation of 
the normalization condition (|37|). Indeed, using the evolution equation for p and the property 
(JBTIj) . one can find 

jt Tlt{p{t)) = (^r ~ ^ ijUij ) Tlt{p{t)) = • (41) 

This property specifies the form of the trace functional up to an overall constant. 

In the Schrodinger picture 3 the physical observables are considered to be chosen once and for 
all at some initial time moment, say t = 0. Only the quantum states evolve according to Eq. 
(I38j) . The expectation value of an observable F(x) at the time moment t relative to a state p(t) 
is given by 

(Fy p = Tr t (F(x)* t p(t,x)), (42) 

where the Wigner function p obeys Eqs. lJBTj) . 

For quadratic Hamiltonians the equation (}3*8*j) coincides with the modified Liouville equation 
(I28j) . In that case, the evolution of an arbitrary quantum state p is governed by the linear first- 
order PDO 

^-~x i n ij { X i i p} + ip,H} = Q, 

for which the initial classical equations (0) play the role of characteristics. So, we arrive at the 
following 

Proposition 1. Let x l (t) = T l -{t)xQ + v l (t) be the general solution to the classical equations 
of motion (QJ) with x l (t)\ t=0 = x l , then the evolution of a quantum state p(t,x) is given by the 
expression 

p(t,x) = p (A(t)[x - v(t)]) , 
where A(t) = r _1 (t) and the initial state po(x) = p(0,x) satisfies Eqs.^Ffy at t = 0. 
3 Of course, all the constructions can be straightforwardly reformulated in the Hciscnbcrg picture. 
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Notice that the quantum evolution of linear systems is completely determined by the classi- 
cal one and does not depend on any ambiguities concerning the choice of the quadratic action 
functional (J2J). 

As we mentioned in Introduction the characteristic feature of dissipation is the presence of 
attractors, i.e. invariant subsets A G M 2n to which all nearby trajectories converge. In other words, 
any classical state supported at a sufficiently small vicinity of A evolves to a state supported at A 
when t — > oo. As the next proposition shows the similar phenomenon takes place at the quantum 
level as well (at least for linear systems with one-point attractor). 

Proposition 2. Given a system of differential equations x = A(t)x having the origin x = as 
the global attractor, then each Wigner's function p(t,x) defines a 5-shaped sequence 



Proof. Since x = is the global attractor lim t _ >00 x(t) = for any solution x(t). On the other 
hand, x(t) = T(t)xo, and hence lim^oo T(t) = 0. According to Proposition 1, for any compactly 
supported function F(x) G C^°(lR 2n ) we have 

J F(x)Mx)d 2n x = -0- J F(x)p (A(t)x)d 2n x = J F(T(t)x)p (x)d 2n x. (43) 

Taking limit, we finally get 

(F>~ = lim J F(T(t)x)p (x)d 2n x = F(0)Tr (p ) = F(0) . 

Thus, 4>t(x) is a 5-shaped sequence. 



4 Examples 

By way of illustration let us consider the deformation quantization of two interesting dissipative 
systems: the damped linear oscillator and the radiating point charge in a homogeneous magnetic 
field. The corresponding classical dynamics has been discussed in Sec. 2, including the time- 
dependent Poisson brackets. In both cases, dissipation has the form of a perturbation over a 
Hamiltonian system and we choose the basis quantum states as the eigen-states for the energy 
(and angular momentum) of the corresponding non-perturbed system. Using these states we then 
consider the evolution for the mean values of energy (and angular momentum) in the presence of 
dissipation. 

4.1 Damped linear oscillator 

Let p(x,p) be the Wigner function describing a pure state of the harmonic oscillator (a = 0) 
with a definite value of energy E. This amounts to saying that p solves the following eigen-value 
problem: 

H * p = p* H = Ep , 

(44) 

P*P = P, P = P, Tr(p) = 1 , 
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where 

H=^(p 2 + u 2 x 2 ) , 

and the *-product (JHHj) is defined at t = by the canonical Poisson brackets {p, x} = 1. Let us 
find p from f)44j) following [27|. The first (complex) equation in ()44|) is equivalent to the pair of 
the real ones 

Hp-£(^ + ?£)=2Ep, {H,p} = 0. (45) 



4 \ dp 2 dx 2 
The second equation implies that 

p = p(H). (46) 

Then the first equation yields 

Hp-^(p"H + p') = -p. (47) 
4 uj 

Introducing new variables 

y = ^H, P = e-y/ 2 f, 

one can bring Eq. (J47j) to the form 

rv + v-v)f f = o. 

This equation is satisfied by Laguerre's polynomials L n (y), provided 



(4* 



E = hjj\n+^\, n = 0,1,2,..., (49) 

and these are known to exhaust all its solutions resulting in integrable Wigner's functions p G 
L 1 (M 2 ). Thus, any eigen- value (}4T?j) corresponds to the unique Wigner's function 

p n {H) = C n exp L n (50) 

The constant C n = (— l) n 2/(j is determined from the normalization condition (|44|) . In Ref. |27j . 
it was shown that the sequence {p n } defines a complete set of orthogonal projectors: 

CO 

Pn*Pm = 5 mn p n , 5(X - x')S(jp ~p') = 1 ^2p n ( y X,p)p n (x',p') . (51) 

rc=0 

Now let us return to the damped linear oscillator with action ([13)1 . The time evolution of the 
mean energy H can easily be calculated using Proposition 1 and the formula (|42J) . We have 

(H)l n = j H(x) Pn (A(t)x)d 2n x = — j H(T(t)x) Pn (x)d 2n x } (52) 

where the initial state p n is one of the states (JSTij) and T(t) solves Eq. (fT2"j) with T(0) = 1. Since 

H(T(t)0 = e- at H(0 , e = (M*. (53) 

we get 

= ^ / dpdxif * p n = " / dpcbp n = e _a *£ n . (54) 

As is seen the quantum evolution of the mean energy relative to the eigen-states p n coincides 
exactly with the classical one: in both cases the energy decreases by exponential law. 
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4.2 Particle in homogeneous magnetic field 

Consider the deformation quantization of the system ()17j) . It is convenient to identify the complete 
set of observables with the energy H and the angular momentum L of the particle without friction 
(A = 0). Then the corresponding set of states with definite values of H and L is determined by 
the equations 

H*p = p*H = Ep, L*p = p*L = Mp, 

(55) 

P*P = P, P = P, Tr(p) = 1 , 



If B \ 2 If B N " 



where E,M G M and 

H = ^[P~ ~^y) + § U + f*J ' L = py-qx. (56) 

The *-product is defined by the formula at t = w.r.t. the canonical Poisson brackets: 
{p, x} = {q,y} = 1 and the other brackets vanish. 

In order to solve (J55j) we introduce the following linear canonical transformation (p, x; q, y) — > 
(P,X;Q,Y): 

B ^ 1 ( B 

P-^y, x =bV + T 

(57) 

1 f B \ B 

together with the functions 

H 1 = P 2 + B 2 X 2 , H 2 = Q 2 + B 2 Y 2 . (58) 

It easy to see that 

H = H U L = B-\H 2 -H 1 ). (59) 

Since Hi and H 2 are nothing but the Hamiltonians of two independent harmonic oscillators, we 
have reduced the eigen- value problem (J55)) to the previous one (I44j) . The Wigner functions solving 
Eqs. flB3|) are given by the (ordinary) products 

Pe,m = Pn{Hi)pi{H 2 ) , (60) 

where p m {Hi,2) is defined by (foTIj) with uj = B, and the eigen- values E and M run trough the sets 

E = hB^n + ^j, M=H(l-n), n, I = 0, 1, 2, ... (61) 

As is seen the eigen-values of the angular momentum L in a state with definite energy E are 
bounded from below by —E/B + h/2. 

Consider now the evolution of the mean values of H and L. Using the fundamental matrix of 
(|19j) and applying Proposition 1 we find 

H(t,t) = H(rm = e 2At H(0, 

(62) 

L(£, t) = L(T(t)0 = + a(t)H(0 + /?(*)#(£) + 7(*M£) , 
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where £ = (x,p,y,qy and 

JV^Q-Py, (63, 

_ Vcos(.Bt).B - 2A 2 e A 'sin( J Bt) - 2A5 
7 ^ ~ 5(52 + ^2) ' 

This gives immediately 

W* B ,m = 7^ / * P E ,m(0 = 7^ y ^Pe,m(0 = e 2M E . (64) 

The mean energy of the particle decreases by exponential low just as it behaves in the classical 
theory. 

The invariance of the Wigner functions Pe,m{Q under reversions in PX- and QV-planes sug- 
gests that (K) e,m = (N) e,m — 0, and hence 

(LY EM = Tr (L(t) * P e,m) = M- a{t)E . (65) 

The same arguments show that 

( x)e,m = ( p)e,m = ( vYem = ( q)e,m = • (66) 

So, at each moment of time the measured values of coordinates and momenta relative to the state 
Pe,m are equal to zero. 

Taking successive limits t — > oo and A — > in (j65)l . we find that the limiting value of the 
angular momentum for a small friction \A\ <C 1 is given by M — E/B. Since (3{t) — > 0, j(t) — > 
as t — > oo, the same limiting value appears in the classical theory as well. 
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